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(N I I. INTRODUCTION 

We consider the model of directed polymers defined in terms of an elastic string 4>{t) directed along the r-axes 
O ■ within an interval [0,t] which passes through a random media described by a random potential V(4>, r). The energy 
of a given polymer's trajectory 0(r) is 
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{\[dr^r)] 2 +VWt),t]}; (1) 



where the disorder potential V[(f>, r] is Gaussian distributed with a zero mean U(</>, r) = and the ^-correlations 
V((j), t)V((J>', t') = uS(t — t')5((J) — 4>') The parameter u describes the strength of the disorder. 

The system of such type as well as the equivalent problem of the KPZ-equation [lj describing the growth in time 
of an interface in the presence of noise have been the subject of intense investigations during the past almost three 
decades (see e.g. |2Hl3{). The breakthrough in these studies took place in 2010 when the exact solution for the free 
energy probability distribution function (PDF) for the model with fixed boundary condition has been found J14l — [21j | - 
It was show that this PDF is given by the Tracy- Widom (TW) distribution of the largest eigenvalue of the Gaussian 
Unitary Ensamblc (GUE) [22]. Since that time important progress in understanding of the statistical properties of the 
t— ( \ KPZ-class systems has been achieved (for the review see J23l.l24|). In particular, by this time it is shown that the free 
energy PDF of the directed polymer model (p} with free boundary conditions is given by the Gaussian Orthogonal 
Ensemble (GOE) TW distribution 0,111], whi le in the presence of a "wall" (4>{t) > 0, < r < t) such PDF is given 
by the Gaussian Simplectic Ensemble (GSE) TW distribution J27|. Besides, the two-point free energy distribution 
function which describes joint statistics of the free energies of the directed polymers coming to two different endpoints 
has been derived in |28l| and quite recently the explicit expressin for the PDF for the end-point 4>(t) fluctuations has 
been obtained [29 32]. 

All these studies, however, describe the statistics of the model (QJ in the so called "one-time" situation. In this 
paper I am going to derive the joint probability distribution function, W(/i, /2, A), for the free energies /i and $2 of 
two directed polymers with fixed boundary conditions at two different times, <fi(ti) — and 4>(t2) — 0, with t\ = i, 
t% = (1 + A)i in the limit t — > oo when the parameter A > remains finite. The derivation is done in terms of the 
Bethe ansatz replica technique. The main points of this approach are described in Section II (for the details of the 
method see e.g. [2(j, [26|, |32|]). Detailed derivation of W(/i, $2, A) is given in Section III with the explicit expression 
of the obtained result presented in Section III.E. Unfortunately, at present stage this result looks rather cumbersom 
and its analytic properties are not clear, although the simplest limit cases can be easily verified, namely: (1) in the 
limit /i — > — oo one recovers the GUE TW distribution for J2\ (2) in the limit fa — > — oo one recovers the GUE TW 
distribution for /i; and (3) in the limit A — > oo one obtains two independent GUE TW distribution for f\ and J2- 
Besides, it has to be noted that the result itself is expressed in terms of specially constructed "generalized" Fredholm 
determinant whose analytic properties, as far as I know, were not yet investigated. All these issues as well as future 
perspectives are briefly discussed in Section IV. 



II. BETHE ANSATZ REPLICA TECHNIQUE 



For the fixed boundary conditions, 0(0) = <f>(t) = x, the partition function of the model ([TJ is 
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exp(-/3F t (x)) 
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where /3 is the inverse temperature and F t (x) is the free energy. In the limit t — > oo the free energy scales as 
(3F t (x) = f3fot + (3x 2 /2t + At/, where /o is the selfaveraging free energy density, At oc t 1 / 3 and / is a random quantity. 
It is the statistics of / which in the limit t — > oo is expected to be described by a non-trivial universal distribution 
W(f). In fact the first two trivial terms of this free energy can be easily eliminated by simple redefinition of the 
partition function, and therefore, to simplify the formulas, these two terms in what follows will be just omitted. 
The calculation of the probability distribution function 



W(x) = lim Prob|7(x) > f] 
is performed in terms of the generating function 
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where (...) denotes the averaging over the random potential V{<j>, r). Instead of the one-point replica partition function 



[Z(x)] one can introduce more general ./V-point object: 



JV 



^(xi,...,XN',t) = Z(xi)Z(x 2 ) .-Z(x N ) = Yl 



<t>a(t) = 



4>a(0) = 



VMr) 



JV 



exp[-/3^ff[0 a ]] 



(5) 



It can be easily shown that \l/(x;£) is the wave function of ./V-particle boson system with attractive ^-interaction: 
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(where k — f3 3 u) with the initial condition ^(x; 0) = H^ =1 5(x a ). The wave function ^(x; t) of this quantum problem 
can be represented in terms of the linear combination of the corresponding eigenfunctions x P q (x) of eq.([6]): 
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Vq |C(q)| 2 * q (x)**(0) exp(-i£(q)) 
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where the symbol J Vq denotes the integration over the momenta (<Zi, •••, #jv)j |C(q)| 2 is the normalization factor and 
E(q) is the eigenvalue (energy) of the eigenstate 'l'q(x). A generic eigenstate of such system is characterized by N 
momenta {q a } (a = 1, ...,N) which split into M (1 < M < N) "clusters" each described by continuous real momenta 
q a (a = 1, ..., M) and characterized by n a discrete imaginary "components" (for details see |33l - l37| ): 
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with the global constraint J2 a =i n a = N. Explicitly, 
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where the summation goes over N\ permutations V of N momenta q a , eq.©, over N particles x a . Note that according 
to the above definition, 'l'q(O) = Nl. The normalization factor corresponding to this eigenfunction is 
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and the eigenvalue: 

N ., M ^ 2 M 
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In this way for the distribution function (0J) one gets the following expression 

(°° (-1") N 1 

i + ESvr exp ( AtAr/ )* (x;i) L a „) (12) 

where the wave function ^(x;i) is defined by eas.(IT))- (|ll[) . Performing the summations over all integers n a ,M and 
N, as well as integrating over all ga in the limit t — > oo one obtains the GUE Tracy- Widom distribution in the form 
of the Fredholm determinant |19l - l2l| : 

W(f) = F 2 (-2- 2 / 3 /) s dct [1 - K_ 2 - 2/3f ] (13) 

where K s is the integral operator on [0, +00) with the Airy kernel: 

K s (uj 1 uj')=j dy Ai(y + u + s) Ai(y + J + s) (w,w'>0) (14) 
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III. TWO-TIME PROBABILITY DISTRIBUTION FUNCTION 

A. Replica Bethe ansatz definition 

Let us consider the situation when one polymer trajectory is arriving to zero at time t\ = t and having (random) 
free energy f tl (0), while the other trajectory is arriving to zero at time £2 = (1 + A)£ and having (random) free energy 
/i 2 (0). One would like to compute the joint probability distribution function 

W(/i,/ 2> A) = lim Prob[/ tl (0) > /1; / t2 (0) > f 2 ] (15) 

In terms of the replica partition functions this quantity can be defined as follows: 
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is the partition function of the directed polymer system in which time goes backwards, from t 2 to t\. For technical 
reasons (for proper regularization of the integration over x at ± infinities) it is convenient to split the partition 
function Z t2 into two parts, the "left" and the "right" ones: 
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Substituting eqs. (ITTl) - (j2"2")) (with t\=t and (t 2 — 1\) = At) into eq. fTo]) and taking into account the definition (J5J) we 
get: 
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where the second (conjugate) wavefunction represent the "backward" propagation from the time moment t 2 — (1 + A)t 
to the previous time moment t\ = t. Schematically the above expression is represented in Figure 1. 




FIG. 1: Schematical representation of the directed polymer paths corresponding to eq. (|23 



B. Calculations 



Substituting the representations (|7|) and © for the wave functions ^(x, 0,y; t) and v I'*(x, y; At) in eq. 
obtain: 
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where to regularize the integrations at ± infinities the supplementary factors ±ie have been introduced (which will 
be set to zero in the final result). The factors Hij(q) {i,j = LK,LR,KR) and Il^ R (p) in the above equation are the 
cross-products of the prefactors of the type [1 — iK/(q^ 
"K", "L" and "R" and "K" and "R" (note that Pkk{<\) 



q(i>)] over the particles belonging to the sectors "L" and 
1). Correspondingly, H* LR {p) is the cross-product of the 



prefactors of the type [1 + iK/(p^> — p™')] over the particles belonging to the sectors "L" and "R". 

The summation over all permutations V of (L + K + R) momenta q over (L + K + R) particles of the sectors "L" 
(xi, ...,xl), "K" (all coordinates equal to zero) and "R" (ya, ...,?/i) can be splitted into four parts: 

(1) the summation over permutations V^ of L (taken at random) momenta q over L particles of the sector "L" ; 

(2) the summation over permutations V^ of R (taken at random) momenta q over R particles of the sector "R"; 

(3) the summation over permutations V^ K ' of K remaining momenta q over K particles of the central sector "K"; 

(4) the summation over permutations , p( LKR ) of the momenta q over sectors "L", "K" and "R". 

Similar splitting can be done for the summation over the permutations V of (L + R) momenta p over (L + R) particles 
of the sectors "L" and "R". In other words, the summations over the permutations in eq. ([23|) can be represented as 
follows: 
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To perform the summations over all these permutations we can use the following important property of the Bethe 
ansatz wave function, cq.®. Namely, it has such structure that for ordered particle's positions in the summation 
over the permutations the momenta components q" belonging to the same cluster also remain ordered. In other 
words, if we consider the momenta of a cluster a, {qf,...,q" }, eq.©, belonging correspondingly to the particles 
{xi t < Xi 2 < ... < Xi na }, the permutation of any two momenta q" and g" of this ordered set gives zero contribution. 
In our case we have three groups of particles: {xi, ..., xl} < {0, ..., 0} < {yn, ..., y\}. Thus, in order to perform the 
summation over the permutations 'P ( ~ LKR ) it is sufficient to split the momenta of each cluster into three parts: 



{q?,...,q:j -> {q?,..;q° n 



( 1„ 



Oi I I Q Q 1 



(27) 



defined by three integer parameters m a , k a and s a , such that m a + k a + s a = n a . In such cluster structure the 
momenta components of the left group {q", ..., q^ } go to the particles of the sector "L"; the momenta components 
of the central group {<7m Q+ i, ■■•,<?„ + k } g° to the particles of the sector "K" ; and the momenta components of the 
right group {<lm +k +1' — > 1m a +k a +s a } S° to the particles of the sector "R". 

Similarly, to perform the summation over the permutations 'p( iH ) of the momenta p each cluster is splitted into 
two parts: {pf, ...,p^ ria \\Pm a +i) ■■■■,Pm a +s a }- Here the components {pf, ...,p^ Q } go to the particle of the sector "L", 
while the components {p" n +1 , ...,p"j +s } go to the particle of the sector "R". In this way, the summations over the 

permutations V^ LKR > and V( LR ) is changed by the summations over the integer parameters {m ai k a , s a } constrained 
by the conditions: 



E 



M 



7 , ka 



K, 



a=l 



M 

E 

a=l 



R 



(28) 



Performing simple integrations over x\, ...,xl and yu, ..., j/i in eq. (|2"5")) and taking into account that all cross-product 



factors Uij are symmetric with respect to the permutations , p( L \'p( K \'p( R \'p( L ) and V^ R \ we get: 
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where Sij is the kronecker symbol and 
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Here the momenta components in the "L" sector are 



q a ->■ q* = q a - -in(m a + k a + S a + 1 - 2r) 



Pa -> p" = -Pa + -ziK(m a + s a + 1 - 2r) 



(31) 
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(where r = 1, ..., m Q ) while in the "R" sector the momenta components are counted "backwards" which makes them 
complex conjugate: 



g« ->■ g*" = g« + -ziK(m a + k a + s a + 1 - 2r) 



p* -► p*" = p« - -i«(m Q + s Q + 1 - 2r) 



(33) 



(where r = 1, ...,s Q ). It turs out that the factors I?i(q, p) and Dfj(q,p), eqs (j30|) -(|3i p . can be nicely represented in 
the determinant form 1381: 



where 
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The expression for D^(q, p) is similar. 

Thus, substituting eq.([29|) into eq. f[24|) and taking into account eqs. ([T0l) - (fTT]) we get: 
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for (m + s) = 0, (m' + s') = 



The structure of the matrix C[p,m + s; p',m' + s'), eq. (|40|) and the factor [(l — <5 m ,o#s,o) + ^(p)^m,o<5s,o] in eq. (|36|) 
take into account the fact that the wave function in the time interval t < r < (1 + A)i exists only provided (m + s) > 1. 
For example, let us consider the cluster structure in which (m a + s a ) > 1 for a — 1, ..., M < M and (m a + s a ) = 
for all a = M + 1, ..., M. In this case for the normalization factor of the wave function VPp(x) (in the time interval 

t < t < (1 + A)t) we get: 
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where M x M matrix C has the structure given in eq.([40|). 
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Now using explicit expressions for the momenta components q" and p", eqs.(|32 |l - (|33|l . we have to express the factors 
Dl, Dr, eqs.(|34 j) - (|35|) . and Q, eq. (|38|l . in terms of analytic functions of the parameters q a ,p a , m a , k a and s a . Due to 
the symmetry of the expressions ([3"o]) and (|3~4"|) - ([3"5]) with respect to permutations of the clusters, the factors Dl and 
Dr in eq.([36f can be represented as follows: 
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Using eqs. ([32"1) - (f33"|) one finds: 



1 



(q? - Pr' - ie ) {i? - p?> - ie + iK ) 

1 



r,r' = l,...,m a 



(c?*™ — p*" + ie) {q*™ — p*" + ie — in) 
d^ = ( iK )" m ° det 



r.r' — l,...,s a 



(42) 
(43) 

(44) 
(45) 



(*«)" 



(x a + r + r' - 2)(x a +r + r' - 1) 

m.!IlSi 2(m °-" 



r,r' — l,...,m 



nS (x Q + Z-l) i ( a ; Q +2m Q -0' 



(46) 



8 



where 
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The expression for D a is obtained by changing m a — > s a and iA a — > — iA a . 
For the factor S^ 1 one obtains the following expression: 
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Using the properties of the Barnes G-function: 



N 



G(N + i) = JJr(i) 



;=i 



G(z + l) = T(z)G(z) 



N 



G{z + N + 1) 



n^+o 0(2+1) 



/=i 

AT 



i=l 



(47) 
(48) 



(49) 



(50) 



(51) 



(52) 



for the factors (|4l)]) and (|49p one eventually obtains the following expressions: 

T{m a + l)G 2 (rn Q + l)G(n a + \k a - m a + iA a ~ ] + 1) G(n a + \k a - m a + iA a ~ } ) 



d a L) = («)- 



G(n Q + ifc Q +iAL } + l) 



G(n a + |fe a — 2m Q + iA a ') 



M 



*? = n 



(iK)- ,n "G(n a + \k a + iA a >)G(n a + \k a - 2m a + iA a ') 
G 2 (m a + l)G 2 {n a + \k a - m a + iA a ~ ] ) 



M 

- n ^ 

a#/3 



where 



(i) _ G(Q a /3 + l)G(Pg/3 + l)G(Q a ^ + m a - mp + l)G(P a/ 3 - m a + mp + 1) 
Q/3 ~ G(Q Q/ 3 + m a + l)G(P a p - m a + l)G(Q a p - mp + l)G(P a p + mp + 1) X 

G{R a p +m a + 171/3 + 2)G(R fja + m a +mp + 2)G(R a p + 2)G(Rp a + 2) 
G(R a p +m a + 2)G{Rp a + mp + 2)G(R a p + mp + 2)G{Rp a + m a + 2) 



(53) 

(54) 

(55) 



Finally, substituting eqs. (|53f and ([54|) into eq. (|42t we get: 
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The expression for Dr is obtained from eq. (|56[) by changing m a — > s a and iA a — > — iA a . In the similar way one 
can derive the analytic expression for the cross-product factor Q, eq.([38|). The final formula for this factor is rather 
cumbersome: it contains the products of all kinds of Gamma functions of the type r[l + i(±m a ± s a ± k a ± mp ± 
sp ± kp) ± i(q a — qp)/ 'k\ and r[l + i(±m Q ± s a ± mp ± sp) ± i(p a — pp)/ k\ (the example of such type of product 
one can see in [32J], eq.(A.17)). We do not reproduce it here as it turns out to be irrelevant in the limit t — > oo (see 
below). 



C. Thermodynamic limit t->oo 

Next step of the calculations is to take the limit t — > oo and to perform the summations over the integers m a , s c 
and k a After rescaling 



la 



2X qa 



p a -> ^ Pa 



with 



A 



l/K 2 i\V3 1 . ., 



5„.2,\l/3 



2V p ) 2 

the expression for the probability distribution function W(fi, /2, A), eq.([55|. can be represented as follows: 

' - (_1)M 

M=l 



W(/l,/ 2 ,A) = lim 1+ V 
M 

>< n 



Q = l 



M! 

dq a dp a 



{2-Kf 



-fc„+s„>l 



m a -\-k a -\-Scx — 1 



(l - <W,o£s Q ,o) + <5(Pa)<5m Q ,0<5s»,0 



x exp 



x det 



M 



|-Afc Q /i - (1 + A) 1/3 A(to q + s a )f 2 - \{m a + k a + s a )ql + -A 3 (m Q + k a + s a ) 

AA(m Q + s a )p 2 a + -A\ 3 (m a + s a ) 3 \ j(-^- 7 ^-,m a ,s a ,k a 

C(q a , \(m a + k a + s a ); qp, X(mp + kp + sp)) det C[p a , \(m a + s a ); pp, \(mp + sp) 



n [b%bw 



J aP ^ap 



a^P 



Here 



J 



(|,|,m,s,fe) = r(m + i)r(s + 1) x 



G(n + \k - m + ^(q - p)(~) + l)G(n + \k 



2A 



(g- P )(+) + l) 



G ( n + \ k + jxii - p) { ~ ] + l ) G ( n + h k -j\(<i- p) {+) + l ) 

G(n+±k+&(q-p)(-))G(n+±k-&(q-p)M) 
G ( n + ifc _ m + £.( q _ p )(-))G(n + §* - 5 - £(« - p)M) 



(57) 



(58) 



(59) 



(60) 
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where n = m + s and (q—p)^' = q — p± ie. 

The summations over the integers m ai s a and k a in the limit A — > oo is performed according to the following 
heuristic algorithm (see also [2g, [32|). Let us consider the example of the series of a general type: 



M 



R= lim n 



a=l 



E (- 1 )™"" 1 



Tl Q =0 



<fr(A; An; n) 



(61) 



where $(A; Azi, ..., Azm; 2i,...,^jy) as a function of the variables {zi, ..., zm} has "good" analytic properties in the 
complex half-plane, Ke(z a ) > 0. Then the summations in the above equation can be represented in terms of the 
following contour integrals: 



M 



R= lim TT 






1 



2i J c sin(7rz Q ) 



$(A; Az; z) 



where the integration contour C is shown in Figure 2a. After rescaling, z a — » z Q /A we get: 

M 



R=U 



1 f dz a 



2ni 



C ^a 



lim $(A; z; z/A) 



(62) 



(63) 



where the parameters z a remain finite in the limit A — > oo. Note that if the summation in eq. (|6ip would start at 
n = 1 (and not at n = 0), the integration in eg. (1631) would go along the contour C shown in Figure 2b. 



-€: 



(a) 











c • 


■••-... 






*\ 





1 2 

(b) 


3 / 



FIG. 2: (a) Contour of integration C in eq. (|62[) ; (b) Contour of integration C' in the case the summation of the series (|6T 
starts at n — 1 



Applying the above algorithm for the expression in eq. (|59p and taking into account that 

lim T(z + 1) = 1 
| 2 j->o 



lim G(z + 1) = 1 

lzl-s-0 



(64) 



we find: 



and 



G ( 2 )||z|«i - z 



lim B 



A^oo 



(L,R) 



lim Q = 1 

A— >-oo 



(65) 



A 11 ? 1 J (!' 

A— >-oo \ A 



P £l £2 £3 
A' A' A' A 



-21 



-22 



Z2 + \zz + ±i(g - p)(-> / 1 z 2 + ±z 3 - ii(g - p)(+) 



J*(q,p, zi,z 2 ,z 3 ) (66) 
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Thus, instead of eg. (|59|) . in the limit A — > oo we get the following much more simple expression: 

M 



W(frJ 2 ,A) = 1+E 



(-1) 



M M 

-n 



M! 

M=l Q=l 



/ + OG /" + 00 7 

d yi ki(y 1+q l-h) Trf^Aif 

-oo J —oo ^* ^ 



y 2 + Ap2_(i + A )i/3 /2 + /l 



AV3 



X exp{(zi Q +Z 2a )(yi +y2) + ^3a2/l} ^(^cnPa^la^a^aJ 



(67) 



x det 



C[q a , (z la + z 2a + z 3a ); qp, (zip + z 2p + z 3 p)) det C[p a , (z la + z 2a ); pp, {zip + z 2 p)) 



where we have introduced the notation 



dq f +0 ° dp 



/ T>(zi,z 2 ,z 3 ,q,p) = — dzx dz 2 — — 8{z 1 )§{z 2 )8{jp) + 

«/ t/ O J « O «* C> *» — OG «/ — oo 

+ [ — [ dzi [ dzJ— - 6(zi)6(z 2 )) 

JC z 3 Jc JC Vz l 2: 2 / 

and used the Airy function identity 

)|-2 3 | = / dy Ai(y) expjzy} 



oo 2tt 7^^ 2tt 



expi 



(68) 



(69) 



D. Generalized Fredholm Determinant 

One can easily see that the expression for the function W(fi, f 2 , A), eq.([67|). has the Fredholm-like structure, 
although it is not the Fredholm determinant. Thus, to proceed further in our calculations we have to introduce a new 
Fredholm-like object. 

In very simplified form the structure we have got in eg. (1671) can be represented as follows: 



M 






M\ 

M=\ a=l 



+ oo 



dx a dy a f(x a ,y a ) 



det 



Ci(x a ,xp) 



det 



a/3=l,...,M 



C 2 (x a ,Xp) 



a/3=l,...,M 



(70) 



Let us introduce the two-dimensional variable £ = (x, y) and define: 



D(Z;S')=D[(x,y);(x',y')] = f(x,y)C 1 (x,a/)C 2 (y,y') 
Then the expression in eg. ((70)) reduces to 



00 / i\M M 



M! 

M=l a=l 



/ 



d 2 ? Q 



det 



D(Z a ;Sp) 



aP=l,...,M 



where 



M 



det 



D(t a ;S p ) 



a0=l,...,M 



EE^ 1 )^- 1 ) 1 ^^'^);^'^)] 



(71) 



(72) 



(73) 



where [V] and [T 7 ] denote the parities of the permutations V and V correspondingly. Thus the object T, eg. ([72)1. can 
be called the generalization of the Fredholm determinant for the two-dimensional variables situation: 



T = det 1 - D 



(74) 
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where D is the integral operator with the kernel given in eq. ([7T|) . One can easily see that just like in the case of the 
usual Frcdholm determinant, the object T can be represented in the exponential form: 



T = expj- ^2 Tr/) 

M=l 



M 



(75) 



where 



TyD 



M 



M 

n 



+oo 



dx a dy a 



D[(xi,yi);(x 2 ,y2)] D[(x 2 ,y 2 y, (0:3,2/3)] D[{x m ,Vm)\ (xi,yi)] (76) 

Coming back to the expression (|57|) , one finds that in terms of the above formalism it reduces to 

00 
W(f x ,h,A) = exp{-]T TrA M } 

where 



(77) 



M=l 



TrA 



M 



M 

n 

Q = l 

e + OO 



P( z 1qi ^2q, z 3a: q a ,p a ) 



-00 J — 00 £i 



i ?/2 + A^-(l + A) 1 / 3 / 2 + / 1 - |- 

" A^ ! 



x exp{(zi Q + z 2a )(yi + y 2 ) + z 3a yi} J*(q a ,p a ,zi a ,z 2a ,z 3a ) 



(78) 



M r 

]J ^ C \9a, (zia + z 2a + z 3a ); q a+l , (z la+1 



+ z 2a+ i + z 3a+1 ) C[p a , (zi a + z 2a ); p a+x , (zi a+ i + z 2a+1 ) 



where, by definition, in the last product it is assumed that the superscript (M + 1) = 1. 

Substituting here eq. (|68[) and taking into account the structure of the matrix C, eq.(|40|). one can easily prove that 
the contributions coming from the first and the second terms of T>(zx a , z 2a , z 3a , q a ,Pa) , eq. (|55jl . completely decouple 
from each other in the trace expression ([75]): 



TrA M = Tr/C 



M 



Tr£ 



M 



where 

Tr/C M = H 
and 



M 



ct=l 



dz. 



3a 



C z 3 



+ OO 



a. J —00 



+00 



-^- J dy Ai(y + ql-fi) exp{yz 3Q } 



M 

n 



1 



z 3a - iq a + z 3a+1 + iq a+ i 



(79) 



(80) 



Tri" = f d -^fdz la fdzj^- -S(z la )S(z 2a) ) r^l r d J>* 
Jc z 3a Jc Jc Kz la z 2a j 



2tt 



2tt 



/+00 1— 

dy 1 Ai(y 1 +ql-f 1 ) / 
CO J —: 



+ °° dy 2 A ./ 2/2 + Ap2_(i + A ) 1 /3 /2 + /l 



AV3 



Ai 



AV3 



x e-xv{{z la + z 2a )(yi + y 2 ) + z 3a yi) J*(q a ,p a , z la , z 2a , z 3a j 



(81) 



M 

n 

a=l 



( z ia + Z 2 Q + z 3a — iq a + z\ a+ i + z 2a+ i + z 3a+1 + iq a +i) 
1 



(zi a + Z 2a + ip a + Z\ a + 1 + Z 2a+1 - ip a+ l) 
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One can easily note that the first contribution, eq. (f80|) . provides the usual GUE Tracy- Widom distribution for the 
free energy A of the directed polymers with fixed (zero) boundary conditions |19T[2l| : 



M 



Tr/C M = Yl 



duj n 



M 

n 



K(u a - A/2 2 / 3 , u a+1 - A/2 2 / 3 ) 



K(uj,uj') = / dy Ai(y + u) Ai(y + u') 



where 



is the Airy kernel. 

Let us consider the second contribution, eq. (|81[) . Substituting 



(82) 



(83) 



and 



(zia + Z 2a + z 3a - iq a + z la+1 + z 2a+1 + z 3a+1 + iq a+ i) = 

= / du} a exp|-(zi Q + z 2a + z 3a - iq a + zi a+1 + z 2a+1 + z 3a+1 + iq a+ i)uj a j 

(z la + z 2a + ip a + z la+1 + z 2a+1 - ip a+1 ) = 

duj a exp|-(zi Q + z 2a - iq a + z la+l + z 2a+1 + iq a+ i)uj a j 



(84) 



(85) 



into eq.([5Tj) we obtain 

M 

Tr£ M = [] 



duj n duj n 



+00 



2tt 



00 j —00 



M ( r-+oo 1 r-+oo 

a 111/ o~ / dyi Ai(yi +g 2 +w Q +w Q _i - A) exp{zg(w Q - w Q _i)} 

V2 t ./ y 2 + Ap 2 +^ + €>«_! - (1 + A) 1 / 3 A + A \ ; . ,- - u 
TJ^ A \ ^TTs J exp{-?p(u; Q -a; Q _i)) x 



x S(q,p,y 1 ,y 2 ) 



(86) 



where 



S 



(q,P,yi,y2) = / — - dz 1 dz 2 ( S(zi)S(z 2 )j exp{(zi + z 2 )(yi + y 2 ) + z 3 yi} 



x 1 



-21 



"« 



zi + ^ 3 -|(g-p)(+) 



Simple integrations over zi, Z2 and 2:3 yield: 

S{q,p,yi,y 2 ) = 47r8(q-p)5{y 1 +y 2 )8(-y 1 ) - %i) % 2 ) - 4% x + y 2 ) %i) 



(87) 



3in[(g-p)yi] 



Substituting eq. 



into eq.dHU) we find 

M 



Tt£ 



A I 



n 



E. The result 



M 



Yl \E[(u> a ,U! a ); (w a+ l,W a+ l)] 



(89) 



where the kernel E(u;u') (with u) = (ui,ui) ) can be represented as follows: 
£[(w,w);(u/,w')] = 
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+ OO 



dq 



dy J ^Ai(y + q 2 + u + u'-f 1 ) Ai 



AV3 

x exp{— iq(ui + uj 1 — uj' — uj) } 



-y + Ag 2 +^ + ^-(l + A) 1/3 / 2 + A 
AV3 



(90) 



A 1 / 3 Jo 
sm[(q-p)y] 

(q-p) 



* r s r°i^+ •*+»+«'- a) *•(■ 



-j/ + Ap 2 +^ + ^-(l + A) 1/3 / 2 + /i 
AV3 



exp{— igfw — u/) + ip(u) — w'J } 



(91) 



1 I™ i [°°j f + °° d( l f + °° d P *■/ 2 , .\ ../y2 + Ap 2 +d; + ^-(l + A) 1 /3/ 2 + / 1 . 



AV3 



x exp|— iq(uj — u/) + ip(u> — uj') } 
Using the Airy function relations: 

I- + QO 



(92) 



A 1 / 3 



dyAi(y + u) Ai(-A 1/3 y + «) = (1 + A)1/3 A 



A 1 ^ 



(1 + A)V3 L(1 + A)V3 

/■ +1 ° ^i Ai( g 2 + M + t-) exp{zg( W -t>)} = 2^ 3 Ai[2^ 3 u] Ai[2 1 / 3 u] 



(93) 

(94) 



redefining (w,w) — >• 2 1 / 3 (1 + A) 1 / 3 (w,a)), and taking into account the contribution (|52"|) . the result for the two-time 
probability distribution function Wyfi, A, A) can be represented in terms of the generalized Fredholm determinant, 
eas.([72" ]) -(|75 ]) . as follows: 



W(hJ 2 ,A) = F 2 (-f 1 /2 2 ^)det[l-£] 



(95) 



Here F2(x) is the GUE Tracy- Widom distribution and £ is the integral operator on [0, +oo) with the kernel 
2?[(w,w); (w',w')l consisting of four contributions: 



£[(«,«); (w'.u/)] - £ 25, [(«,«); («',«')] 



(=i 



where 



£i[(w,w);(w',w')] = AiL + w'-/ 2 /2 2 / 3 ] AiL'+w-/ 2 /2 2 / 3 ] 



(96) 



(97) 



l + A\2/3 



E 2 [( U ,Q);(u;',Q')} = -2 1/3 (^) 






+ A- 2 /V + (1 + A) 1 / 3 ^ W) - A 



x Ai 



-y + A 1 / 3 ? 2 + (1 + A)V3(£ + £') - (1 + A)V3/ 2 + / a 



AV3 



•1 + A\ i/3 



x exp< — iq\ — - — I [u! + uj — uj — us) > 



(98) 
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1 + A\2/3 



2% [(«,«); (c/.fi/)] = -^ /3 {^) 



dy 



+ °° dqdp .. 
Ai 



x Ai 



J J-oo ( 27r 

y + A^V + (1 + A) 1 /3( ci + £') _ (i + A)V3/ 2 + ^ 



jy + A- 2 / 3 g 2 + (1 + A) 1 / 3 ^ + «') - A 



A 1 / 3 



sin [A 1 / 3 (g — p)y] 
(<7-p) 



cxp 



'1 + A\V3 



/l + A\l/3 / 






(99) 



£ 4 [( W ,w);( W ',£')] 



(1 + A) 2 / 3 
A 1 / 3 



dyi dy 2 Ai 



2/i 



+ (1 + A) 1/3 cj - A/2 2/3 



Ai 



x Ai 



1 + An V3. (l + A) 1 / 3 ^-/! 
2/2 + I — : — ) w 



22/3^1/3 



Ai 



2/! + (1 + A) : /V - A/2 2 / 3 x 

(100) 



1 + Ax i/3 w (l + A) 1 ^-/! 
2/2 + I — : — ) uj 



22/3^1/3 



IV. CONCLUSIONS 



The result for the two-time free energy distribution function, eq. ([95|) - (ll00|) , obtained in this paper should be 
considered as the prelimenary one. At present stage it is expressed in terms of specially constructed generalized 
Fredholm determinant, eqs.([72 |) - ([76l) . whose analytic properties are not yet studied. Besides, the the final expression 
for the kernel of the integral operator i?(a;;u/) defining the distribution function, eq. (|9"5|) . is rather cumbersome and 
the study of its analytic properties requires special investigation which will be done elsewhere. At present stage only 
three limit cases for the distribution function W(f\, f 2 , A) can be easily checked: 

(1) The limit /1 — > —00. In this limit one easily finds that F 2 [— /i/2 2 / 3 ] — > 1, and the three contributions 
-^2,3,4 (<■»■';<■»■>')' ec l s -(EHl)- (HQUJ) , turn to zero, so that the only contribution which survive is -Ei(u>;u>') given in (|97p . 
Substituting its expression into the integral representation of the Fredholm determinant exp{— X)m=i Tr£ M f one 
easily recovers the usual Fredholm determinant with the Airy kernel. In other words, 



lim W(/i,/ 2 ,A) = F 2 [-f 2 /2 2 ^ 3 ] 



(101) 



(2) The limit f 2 — > —00. In this case all four contributions (|97l) - (|100j) turn to zero, so that 

lim W(hJ 2 ,A) = F 2 [-h/2 2 ' 3 ] 

/ 2 ->-00 



(102) 



(3) The limit A — > 00. In this case the three contributions (|98|) - (|100p turn to zero, so that one finds: 

lirn^ W(h,f 2 ,A) -.-. F 2 [-f 1 /2 2 / 3 ]F 2 [-f 2 /2 2 ^ 3 ] ill):]) 

as it should be. 



A->oo 



The limit A — ¥ is much more tricky. First of all, technically it is not so easy to study, and second it is not quite 
clear what kind of behavior for the probability distribution function should be expected in this case, as the result 
(l95l) - (|100[) gas been derived in the limit (£2 — ii) = Ai — >• 00. In this situation the physical meaning of the limit 
A — > is not evident. 

Note finally that the obtained result can be easily generalized for the case when the directed polymer at time t\ 
comes not to zero but to a given point x\, while at time t 2 it comes to another given point x 2 . After the proper rescaling 

x \,i - ► (f3ut 2 ) Xi t 2 the expression for the corresponding distribution function W(fi, Xi, f 2 , x 2 , AJ is obtained from 
W(fi,f 2 , A) by the trivial shift: 



W(h,x x ,f 2 ,x 2 ,A) = W(f 1 +x 2 /2,f 2 + x 2 /2,A) 



(104) 
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